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It is shown that under changing density or temperature a nucleon Fermi superfluid can undergo 
a phase transition to an anisotropic superfluid state, characterized by nonvanishing gaps in pairing 
channels with singlet-singlet (SS) and triplet-singlet (TS) pairing of nucleons (in spin and isospin 
spaces). In the SS pairing channel nucleons are paired with nonzero orbital angular momentum. 
Such two-gap states can arise as a result of branching from the one-gap solution of the self-consistent 
equations, describing SS or TS pairing of nucleons, that depends on the relationship between SS 
and TS coupling constants at the branching point. The density/temperature dependence of the 
order parameters and the critical temperature for transition to the anisotropic two-gap state are 
determined in a model with the SkP effective interaction. It is shown that the anisotropic SS-TS 
superfluid phase corresponds to a metastable state in nuclear matter. 
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I. INTRODUCTION 

It is known that at sufficiently low temperatures a nu- 
cleon Fermi system becomes unstable with respect to for- 
mation of Cooper pairs due to the attractive component 
of the nucleon-nucleon (NN) potential. Basically in early 
articles on superfluidity isovector pairing of the same nu- 
cleons in nuclei and nuclear matter ^ was studied. 
Later it was realized that isoscalar neutron-proton pair- 
ing plays an essential role in the description of superflu- 
idity of finite nuclei with N ~ Z (see Refs. (3[ || and 
references therein) and nearly symmetric nuclear mat- 
ter ||— Possible coexistence of isoscalar and isovector 
nucleon pairing in finite nuclei was studied in Ref. 0. 
In addition to isovector pairing, isoscalar pairing modes 
have been investigated to explain the excitation ener- 
gies in N — Z nuclei |t0|. As was shown in Ref. [ pd| , 
at low densities of nuclear matter the coupling between 
isospin T — and T = 1 pairing channels may be of 
importance, leading to the emergence of multi-gap su- 
perfluid states with nonvanishing gaps in both pairing 
channels. In Refs. jyj |§ the case of a multi-gap con- 
densate was considered, when the energy gaps in 1 Sq and 
3 Si pairing channels are nonzero. For such a conden- 
sate the orbital angular momentum L of a pair is equal 
to zero and, hence, the order parameters in both pair- 
ing channels arc isotropic functions of the momentum. 
However, we can consider another possibility, when the 
pairing of nucleons in one (or both) pairing channels oc- 
curs in a state with nonzero orbital angular momentum. 
In this case an anisotropic multi-gap condensate will be 
realized, where the order parameter(s) will depend not 
only on the absolute value of momentum, but also on its 
direction. Further we will consider an isoscalar multi- 
gap superfluid state, corresponding to a superposition 
of pair states with spin 5 = and S = 1. According 



to the Pauli principle, the orbital angular momentum in 
the first pairing channel (S = 0, T = 0) is odd and in 
the second one (S = 1,T = 0) it is even. Hence, such 
a multi-gap condensate will necessarily be anisotropic. 
Our main purpose will be the clarification of a mecha- 
nism for the appearance of anisotropic multi-gap super- 
fluid states, finding the corresponding critical temper- 
ature and comparing the free energies for different su- 
perfluid phases. As a theoretical framework, we use the 
Fermi-liquid (FL) approach 13 , in which normal and 
anomalous FL interaction amplitudes are treated on an 
equal footing. As a NN interaction we choose the density 
dependent Skyrme effective forces, used earlier in a num- 
ber of contexts for calculations in finite nuclei ju| |l5| and 
infinite nuclear matter ]il| -|T9l . In specific calculations 
we use the SkP potential |14[ , for which the strongest 
interaction in the state of a Cooper pair with nonzero or- 
bital angular momentum is realized in the pairing channel 
S = 0, T = 0, and in the state with L = it is in the 
channel S = 1,T = 0. 



II. BASIC EQUATIONS 

Superfluid states of nuclear matter are described 
by the normal f K1K2 = TV ga^ a Kl and anomalous 
<?kik 2 — Tr ga K . 2 a Kl distribution functions of nucleons 
(k = (p, a, t), p is momentum, <t(t) is the projection 
of spin (isospin) on the third axis, g is the density ma- 
trix of the system). The energy of the system is speci- 
fied as a functional of the distribution functions / and 
9, E = E(f,g). It determines the quasiparticle energy e 
and the matrix order parameter A of the system 
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The self-consistent matrix equation for determining the 
distribution functions / and g follows from the minimum 
condition of the thermodynamic potential [|l3[ and is 



Here 



/ = {cxp(r e + Y 4 ) + l} 1 = {exp(Yo|) + l}~ 
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Here the quantities e, A, I4 are, in turn, matrices in 
the space of the k variables, with 1^ K1K2 — Y4 T1 S K1K2 
(ti = p,n), Y Q = l/T, Y ip = -fip/T and Y in = -<U„/T 
are the Lagrange multipliers, fi p and /i n are the chemical 
potentials for protons and neutrons, T is the tempera- 
ture. We shall study two-gap superfluid states in sym- 
metric nuclear matter, corresponding to a superposition 
of pair states with total spin S and isospin T: S = 0, 
T = (singlet-singlet (SS) pairing in spin and isospin 
spaces) and S = 1,T = (triplet-singlet (TS) pairing) 
with the projections S z = T z = (SS-TS states). In this 
case the normal / and anomalous g distribution functions 
read |§ 

/(p) = /oo(p)ctot , (3) 

giP) = 0Oo(p)c 2 T 2 +330(P)0'30'2T2, 

where <7j and are the Pauli matrices in spin and isospin 
spaces, respectively. The components of the anomalous 
distribution function in Eq. (|^) possess different symme- 
try properties, 



0oo (~p) = -Soo(p), 03o(-p) = 030 (p), 



(4) 



and, hence, considering this, the multi-gap condensate 
will be anisotropic. For the energy functional, which is 
invariant with respect to rotations in spin and isospin 
spaces, the structure of the single particle energy and 
the order parameter is similar to that of the distribution 
functions /, g: 

e(p) = £oo(p)ctoTo, (5) 

A(p) = A o(p)er 2 T2 + A 3 o(p)cr 3 cr2T2, 



where 



Aoo(-p) = -Aoo(p), A 30 (-p) = A 30 (p) 



(6) 



Using the procedure of block diagonalization p3| , 
one can evidently express the distribution functions 
/ooj0ooi03o m terms of the quantities e and A: 
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E ± = VC 2 + |A±| 2 , A±=Aoo±A 30 , 

£(p) = £oo(p) - Mo, n± = {exp(y £±) + 

Ho is the chemical potential, which should be determined 
from the normalization condition 



(10) 



g is density of symmetric nuclear matter. As follows from 
Eqs. the nucleon superfluid is characterized by 

two types of fermion excitations with gaps A± in the 
spectrum. In this case the spectrum is two-fold split 
due to coupling of SS and TS pairing channels (Aoo ^ 

<).A; i:: i 0). 

To obtain the self-consistent equations for the quanti- 
ties A and £, it is necessary to specify the energy func- 
tional of the system, which we write in the form 

E(f,g) = E Q (f)+E mt (f)+E mt (g), (11) 

Eo(f) = 4]>>o(p)/oo(p), E mt (f) = 2^e o(p)/oo(p), 
p p 

£o(p) = £oo(p) = u o(k)foo{q), k = 



2V 



2 t 

Eint{g) = ^Xl(5oo(p) 1/ b(p,q)0oo(q) 

p,q 

+ 03o(p) y i(P:q)03o(q), 



Here mo is the bare mass of a nucleon, f/o(k) is the nor- 
mal FL amplitude, Vb(p, q), Vx(p, q) are the anomalous 
FL amplitudes, describing interactions in the SS and 
TS pairing channels, respectively. With allowance for 
Eqs. ([!]), (fil"|), we obtain self-consistent equations in the 
form 

£(P) =£q(p) -Mo + £oo(p), (12) 
Aoo(p) = y. J2 ^o(P, q)0oo(q), (13) 



A 30 (p) = y53Vi(p,q)03o(q). 



(14) 



Taking into account Eqs. (g), (g), we obtain equations 
for the energy gaps A 00 , A 30 
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Eqs. (|l2l), (|l5|), (|l(jl ) describe two-gap superfiuid states of 
symmetric nuclear matter and contain one-gap solutions 
with A 00 ^ 0, A 30 = (SS pairing) and A 00 = 0, A 30 ^ 
(TS pairing) as some particular cases. 

To obtain numerical results we will use the Skyrme 
effective interaction, for which the normal and anomalous 
FL amplitudes read fl8[ 



U (k) = 6t Q + t 3 g 
+ £[3* 
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where ti,Xi,/3 are phenomenological constants, charac- 
terizing the given parametrization of Skyrme forces. (W e 
shall use the SkP 1^] potential.) According to Eqs. (18), 
(|l9|), pairing in the SS channel occurs with orbital an- 
gular momentum L — 1 and in the TS channel with 
L = 0. Note that in the case of the effective Skyrme 
interaction the normal FL amplitude Uq is quadratic in 
momentum and hence describes a renormalization of free 
nucleon mass and chemical potential. The expression for 
the quantity £, given by Eq. (|l2]), with regard for the 
explicit form of the amplitude Uq, reads 
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where the effective nucleon mass m is defined by the for- 
mula 
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and the effective chemical potential /j, should be deter- 
mined from Eq. (TOT 



III. 



ORDER PARAMETERS AT ZERO 
TEMPERATURE 



Let us consider the case of zero temperature. We shall 
analyze Eqs. (15),(|l6|), using the simplifying assumption, 
that FL amplitudes Vo, 14 are nonzero only in a narrow 
layer near the Fermi surface: |£| < 9, 9 <C ep (we set 
9 = O.lep). Then the TS energy gap represents some 
constant quantity A30 = A 3 o(p ~ Pf), while the SS en- 
ergy gap is angular dependent, Aoo = Ap°n, n being an 
arbitrary real unit vector. In addition, we shall neglect 
the influence of the finite size of the gaps on the chem- 



ical potential fj, and set /i = 



2m 



kt 



1/3 



As 



a result of these assumptions, we arrive at equations for 



determining the quantities A, A30 
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|Ax ± A 30 | 2 , g 0A = -v f Vo A (pf,Pf) 



where 
E± = V? 

and vp — zpink is the density of states at the Fermi 
surface for a nucleon with the given spin and isospin pro- 
jections. Our main goal is to find two-gap solutions with 
A ^ 0, A30 7^ and to clarify the mechanism of their 
appearance. To give some analytical consideration, we 
will assume that the conditions A, A30 <C 9 are fulfilled 
(logarithmic approximation). As will become apparent, 
this approximation works quite well just in the density 
region where the two-gap solutions exist. Introducing 
the ratio of the energy gaps a = A30 / A and performing 
integration in Eqs. (p0|),(21 ) in the logarithmic approxi- 



mation, we arrive at equations for the quantities a and 
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Excluding A from Eqs. (|22|), we obtain the equation 
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Since tp m in = y(0) = — 1/3 (the point a = is the 
point of removable discontinuity) and <p max = y(±oo) = 
1, then Eq. (|2^) has a solution for a, if the coupling 
constants g$ and g\ satisfy the inequalities 
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(24) 



These restrictions have to be fulfilled in the logarith- 
mic approximation for the existence of the SS-TS mixed 
state. The sense of the restrictions (24) on SS and TS 



coupling constants is that the quantities gi and <?o/3 
must be of the same order of magnitude. Clearly, sim- 
ilar restrictions exist in a general case, when the condi- 
tions of the logarithmic approximation are not fulfilled. 
Note that the constraints ( |24| ) are more strict than the 
corresponding ones for the existence of a TS-ST mixed 
state (n| :-l < l/g x - l/g 2 < 1, g 2 = -vfV 2 (pf,Pf) 
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FIG. 1: Order parameters A, A30 vs. density at zero tem- 
perature for SkP force; ss(mix) and is(mix) are notations for 
the dependencies of the energy gaps A and A30 in the mixed 
SS-TS solution. 



(the constants 31,(72, determined here, are four times as 
large as those in Ref. |p"lf). 

The results of a numerical determination of the energy 
gaps A(g), A 30 (g) as functions of density from Eqs. (20), 
( pl| ) are presented in Fig. [j], where the part of the phase 
diagram corresponding to the mixed SS-TS states is 
shown. One can see that SS-TS solutions exist in the 
finite density region g± < g < g 2 . In the left critical 
point two-gap solutions appear as a result of branch- 
ing from a one-gap SS solution (at the branching point 
A 30 (f?i) = 0, A(ei) = A ss ( gi ), A ss being a one-gap SS 
solution). At g — g\ the coupling constants are related 
by the formula 



1 



3i(0i) 9o(Qi 



(25) 



In the logarithmic approximation the one-gap SS and TS 
order parameters have the form 



in the density domain where the two-gap solutions exist: 
the maximum value of the ratios A/9, A 3 o/# does not 
exceed 0.26. 

From Fig. [I] it is seen that SS-TS mixed states exist in 
a density interval that is much closer to nuclear matter 
saturation density than that for TS-ST multi-gap states 
pit , which exist in the region g < g', where g 1 w 0.05 -r- 
0.06 fm~ 3 . Hence, there is no competition between SS- 
TS and TS-ST superfluid states, which exist in quite 
different density domains. 



IV. CRITICAL TEMPERATURE, ORDER 
PARAMETERS AT NONZERO TEMPERATURE 



The analysis given in the previous section relates to 
the case of zero temperature. It is clear that if SS-TS 
states exist at T=0, then such states appear first at some 
critical temperature. To determine the critical temper- 
ature we use the following considerations. Obviously, 
SS-TS solutions arise as a result of branching from SS 
or TS one-gap solutions of Eqs. (^5|),(|l^). If branch- 
ing occurs from a TS solution then at the critical point 
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A -» in Eqs. (|K 
mining T c 
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we obtain equations for deter- 
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The first of these equations determines the temperature 
behavior of TS energy gap, the second one determines 
the critical temperature T c at which the mixed SS-TS 
solution branches from the one-gap TS solution. If SS- 
TS states appear from a SS solution, then A 30 (T C ) = 
0,A(T c ) = A SS (T C ). Considering the limit A 30 -> in 
Eqs. (|l5|),(16), we obtain 



/ 3 In , /lN l = g d£ dx — tanh , E = j£ 2 + A 2 x 2 , 

A ss = 20exp( -- + -), A 3 s = 2^exp( ). y ° J q J E 2T C ' V? 

v 3o 3/ V 31 ' , r 



From here and Eq. ( P5| ) it follows that at zero tempera- 
ture at the left branching point the energy gaps in the SS 
and TS pairing channels are equal, A ss (gi) — A ts (gi). 
This peculiarity is preserved for finite temperatures, as 
will be seen in Section IV. 

In the right critical point the two-gap solutions branch 
off from a one-gap TS solution (at the branching point 
A(g 2 ) = 0, A 30 {g 2 ) = A 3 s (p 2 )). At g = g 2 
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Note that the conditions of the logarithmic approxima- 
tion, A <c 9, A 30 9, are fulfilled quite satisfactorily 



Here the first equation determines the temperature be- 
havior of the SS energy gap, the second one determines 
the critical temperature T c at which the mixed SS-TS 
solution branches from the one-gap SS solution. 

The results of a numerical solution of Eqs. (p7[),(|2"s|) 
and Eqs. (p9|) , (|30|) are presented in Fig. ||. One can see 
that the curve T c (g) consists of two branches. The left 
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FIG. 2: Critical temperature of SS-TS superfluid state vs. 
density for SkP force; notations ss and ts correspond to 
branching of SS-TS solutions from SS and TS one-gap so- 
lutions, respectively. 



branch corresponds to the appearance at a critical tem- 
perature T c of a SS-TS solution from the SS one-gap so- 
lution, the right one corresponds to the appearance of a 
SS-TS solution from the TS one-gap solution. The max- 
imum value of T c is approximately equal to 0.38 MeV at 
density g m ~ 0.139 fm~ 3 . In the limit T c — > 0, from 
Eqs. @,(||), we obtain Eq. @ for the right critical 
point q = q 2 (Q2 ~ 0.15 fnT 3 ) and from Eqs. ©,© 
we obtain Eq. p5| ) for the left critical point g = gi 
(qx rj 0.134 fm'^jTThus, in the density interval (qi, g m ) 
SS-TS solutions appear as a result of a phase transition 
in temperature from a one-gap SS solution, and in the 
interval (g m , Q2) they appear from a one-gap TS solu- 
tion. If gi < g < g m , the coupling constants satisfy the 
inequality g\ > <?o/3, and for g m < g < 02 it is g\ < go/3. 

To determine the temperature behavior of the order 
parameters, one should consider Eqs. (|i~5|), (|l6|). Ac- 
cording to our analysis we can consider two possibili- 
ties, when branching occurs at density g such that (1) 
gi < g < g m and (2) g m < g < g%. The results of nu- 
merical calculations are shown in Fig. ||. In the first case 
(Fig. |(a)) we have A 30 (T C ) = 0, A(T C ) = A SS (T C ), and, 
as in the case of zero temperature, the one-gap order 
parameters in the SS and TS pairing channels are also 
equal, A SS (T C ) = A^ (T C ). In the second case (Fig. |(b)) 
A(T C ) = 0,A 30 (T c ) = A|u(T c ). Thus, the temperature 
region T < T c corresponds to anisotropic multi-gap su- 
perfluidity when, together with one-gap solutions, we 
have two-gap solutions with nonzero SS and TS order 
parameters in both pairing channels. 
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FIG. 3: Order parameters A, A30 vs. temperature at density 
(a) g = 0.138 fm" 3 and fb) g = 0.140 fm" 3 . Other notations 
are the same as in Fig. hi 



V. THERMODYNAMIC STABILITY 

Since we have a few solutions of the self-consistent 
equations it is necessary to check which solution is ther- 
modynamically favorable. For this purpose it is necessary 
to determine the free energy of the corresponding states. 
It consists of two terms, F = E(f,g) — TS(f,g), where 
S is entropy of the system. Taking into account (u])-(||), 
the energy functional (O) is 
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e(p) = ^-,E ± = V^ 2 + |A±| 2 , A± = A 00 ± A 30 . 

The entropy of the system, given in a general theory of 
superfluid FL by the expression S = — Tr/ln/ |13|, can 
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FIG. 4: Free energy density, measured from that of the nor- 
mal state, for SS, TS and SS-TS superfluid states at zero 
temperature. 



be represented in the form 

S = -2 E P (n+ In n+ + (1 - n+) ln(l - 



+n_ In n_ + (1 — n_) ln(l — , 

where n± = {exp(E±/T) + The results of a nu- 

merical calculation of the free energy density, measured 
from that of the normal state, for the case of zero tem- 
perature, are given in Fig. [|. One can see that in the 
density interval g\ < g < g rn the TS superfluid phase 
is thermodynamically most preferable as compared with 
other phases, and at g m < g < gi the SS superfluid 
state wins competition for thermodynamic stability. In 
both cases the mixed SS-TS state appears as a result 
of a phase transition in density from a thermodynam- 
ically less favorable one-gap superfluid state (SS state, 
if gi < g < g m , and TS state, if g m < g < Q2) and 
corresponds to a metastable state in superfluid nuclear 
matter. In Fig. || we show the difference between the 
free energy densities of superfluid and normal states as 
a function of temperature. Fig. |^(a) corresponds to the 
branching of the SS-TS mixed state from the one-gap 
SS solution at fixed density in the range Qi < g < g m , 
and Fig. ||(b) depicts branching from the one-gap TS so- 
lution at a density in the interval g m < g < Qi. As seen, 
for all temperatures T < T c , the SS-TS superfluid phase 
corresponds to a metastable state in superfluid nuclear 
matter. 



CONCLUSION 

We have considered the possibility of the formation of 
an anisotropic multi-gap condensate in superfluid sym- 
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FIG. 5: Difference between free energy densities of super- 
fluid and normal states vs. temperature at density (a) g = 
0.138 fm -3 and (b) g = 0.140 fm -3 . 



metric nuclear matter, corresponding to the superposi- 
tion of states with SS and TS pairing of nucleons. In the 
SS channel, pairing occurs with nonzero orbital angular 
momentum and hence the energy gap is an anisotropic 
function of momentum. The self-consistent equations for 
such two-gap states differ essentially from the equations 
of BCS theory and contain one-gap solutions (SS and TS) 
as particular cases. The analysis of the self-consistent 
equations at zero temperature in the logarithmic ap- 
proximation shows that anisotropic multi-gap superfluid 
states can exist only under quite specific restrictions on 
the coupling constants go and g\, describing interaction 
of nucleons in the SS and TS pairing channels. Since the 
constants of the effective interaction depend on density, 
there are density domains, where one-gap or anisotropic 
two-gap solutions exist. Calculations with the effective 
SkP interaction, chosen as the model of the NN inter- 
action, indicate that two-gap SS-TS states can arise in 
nuclear matter as a result of a phase transition in density 
from a one-gap SS or TS state. In the first case at critical 
density gi > go/3, in the second, the opposite inequal- 



7 



ity is valid. Comparing free energies, branching occurs 
from thermodynamically less favorable one-gap solution 
and hence the anisotropic two-gap supcrfluid state cor- 
responds to a metastable state in nuclear matter. Deter- 
mination of the critical temperature T c of the transition 
to the SS-TS state as a function of density shows that 
the corresponding curve consists of two branches. One of 
them is related to the appearance of a SS-TS anisotropic 
state as a result of branching at T c from the one-gap SS 
solution, another is related to branching from the one- 
gap TS solution. Studying the temperature behavior of 
the order parameters shows that mixed two-gap solutions 
exist for temperatures T <T C . Comparison of free ener- 



gies leads to the conclusion that the anisotropic SS-TS 
phase represents a metastable state for the whole tem- 
perature interval T < T c . Calculations show that mixed 
SS-TS states exist in a density domain, that is close to 
nuclear matter saturation density, in contrast to TS-ST 
mixed states, which only exist in the low density domain 
of nuclear matter. 

Acknowledgement. A.I. is grateful for the hospital- 
ity and support of Rostock University, where a signifi- 
cant part of this work was completed. Discussions with 
S. Peletminsky and A. Yatsenko are gratefully appreci- 
ated. A.I. acknowledges the financial support of STCU 
(grant No. 1480). 



[1] A. Bohr and B.R. Mottelson, Nuclear structure (Ben- 
jamin, New York, 1969), Vol. 1. 
[2] P. Ring and P. Schuck, The Nuclear Many-Body Problem 

(Springer, New York, 1980). 
[3] A. Goodman, Nucl. Phys. A352, 30 (1981); A352, 45 

(1981); A369, 365 (1981). 
[4] G. Ropke, A. Schnell, P. Schuck, and U. Lombardo, Phys. 

Rev. C 61, 024306 (2000). 
[5] Th. Aim, B.L. Friman, G. Ropke, and H. Schulz, Nucl. 

Phys. A551, 45 (1993). 
[6] M. Baldo, U. Lombardo, and P. Schuck, Phys. Rev. C 

52, 975 (1995). 
[7] Th. Aim, G. Ropke, A. Sedrakian, and H. Schulz, Nucl. 

Phys. A604, 491 (1996). 
[8] A. Sedrakian, T. Aim, and U. Lombardo, Phys. Rev. C 

55, R582 (1997). 
[9] A. Goodman, Phys. Rev. C 60, 014311 (1999). 
[10] W. Satula and R. Wyss, Phys. Rev. Lett. 86, 4488 (2001); 

87, 052504 (2001). 
[11] A.I. Akhiezer, A. A. Isayev, S.V. Peletminsky, and A. A. 



Yatsenko, Phys. Lett. B 451, 430 (1999). 
[12] A. A. Isayev, Phys. Rev. C 65, 031302(R) (2002). 
[13] A.I. Akhiezer, V.V. Krasil'nikov, S.V. Peletminsky, and 

A. A. Yatsenko, Phys. Rep. 245, 1 (1994). 
[14] J. Dobaczewski, H. Flocard, and J. Treiner, Nucl. Phys. 

A422, 103 (1984). 
[15] P.-G. Reinhard, D.J. Dean, W. Nazarewicz, J. 

Dobaczewski, J. A. Maruhn, and MR. Strayer, Phys. 

Rev. C 60, 014316 (1999). 
[16] R.K. Su, S.D. Yang, and T.T.S. Kuo, Phys. Rev. C 35, 

1539 (1987). 

[17] M.F. Jiang and T.T.S. Kuo, Nucl. Phys. A481, 294 
(1988). 

[18] A.I. Akhiezer, A.A. Isayev, S.V. Peletminsky, A.P. 

Rekalo, and A.A. Yatsenko, Zh. Eksp. Teor. Fiz. 112, 

3 (1997) [Sov. Phys. JETP 85, 1 (1997)]. 
[19] A.I. Akhiezer, A.A. Isayev, S.V. Peletminsky, and A.A. 

Yatsenko, Phys. Rev. C 63, 021304(R) (2001). 



